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Abstract 

We propose a cosmological model of D3-brane universe on compact and stable supergravity back- 
ground of wrapped D7-branes in type IIB string theory previously argued to be dual to pure J\f = 1 
SU(N) gauge theory in four dimensions. A model universe of order Planck size near the UV bound- 
ary dynamically flows toward the IR with constant total energy density and accelerating expansion 
followed by smooth transition to decelerating expansion and collides with the wrapped D7-branes 
at the IR boundary. The model addresses the horizon and flatness problems with most of the 
expansion produced during the decelerating expansion phase. The inflationary scenario is used to 
generate sources of inhomogeneities in the cosmic microwave background radiation and seeds for 
large scale structure formation from quantum fluctuations which exit the Hubble radius early during 
the accelerating expansion phase and the model addresses the inhomogeneity problem with red tilt 
in the power spectrum. We propose that the kinetic energy of the model universe is converted to 
matter and radiation by the collision followed by formation of baryons that stabilizes the model 
universe against gravitational force from the background at a finite distance from the IR boundary 
with the wrapped D7-branes serving as sources of color. Friedmann evolution then takes over with 
a positive cosmological constant term coming from the remaining potential energy density which is 
interpreted as dark energy. The magnitude of dark energy density is smaller than the total energy 
density during the flow by a ratio of the scale factor when the model universe appears in the UV 
to the scale factor at the moment of collision and stays constant while the matter-radiation density 
falls during Friedmann expansion. 
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1 Introduction 

One important experimental observation in early 20 th century cosmology was that the uni- 
verse was not static but expanding [1], and according to the hot big bang hypothesis the 
universe originated with sub-Planckian size and expanded to its current size. Two of the 
limitations of the old hot big bang model are the horizon problem (on why causally discon- 
nected regions look nearly identical) and flatness problem (on why the universe is spatially 
nearly flat, which requires a highly fine-tuned initial rate of expansion for a given initial 
matter distribution). 

The inflationary paradigm [2-4] addresses these limitations by positing that the uni- 
verse underwent a period of accelerated expansion whereby an initially Planck size universe 
gets exponentially stretched by a factor of about 10 28 or more which was then followed by 
Friedmann evolution of the old big bang model. In standard inflationary models, inflation 
corresponds to the motion of a scalar inflaton field in a given potential, and a suitable inflaton 
potential is introduced by hand and its consequences are investigated which, in addition to 
resolving the limitations of the old big bang model, serve as a probe to fundamental physics 
at high energies through current cosmological observations of the early universe. At the end 
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of inflation, the scalar field is believed to decay to particles which fill the early universe in 
a process called reheating. As the universe expanded and cooled, baryons, atoms, nuclei, 
and then later large scale structures whose origins are quantum fluctuations during inflation 
formed. The same quantum fluctuations are sources of thermal fluctuations in the cosmic 
microwave background radiation (CMBR). The universe is nearly homogeneous and isotropic 
with thermal fluctuations in the CMBR being of order 10~ 5 . Details in the power spectrum 
together with required magnitude of inflation needed to address the horizon and flatness 
problems put constraints on the inflaton potential. 

A second more recent important cosmological observation was that the expansion of the 
universe is currently accelerating [5], which is attributed to the presence of dark energy 
accounting for about 70% of the total energy content of the universe. 

Flux compactifications in string theory provide a scheme for constructing models of the 
universe in which the cosmic evolution could be dynamically determined for a given back- 
ground. A brane inflation model in which the universe is a D3-brane and the inflaton field is 
the position of the D3-brane moving due to interaction with an anti-D3-brane was proposed 
in [6] and further investigation on stabilized supergravity background was initiated in [7]. 

A duality between N D7-branes filling four-dimensional spacetime and wrapping a 4-cycle 
on R 1 ' 3 x || x y~ x ^ and pure Af = 1 SU (N) gauge theory in four dimensions was proposed 
and argued in [8]. The wrapped D7-branes turn on all F\, F3, H3, and F5 fluxes of type IIB 
theory and induce torsion. The supergravity solutions were explicitly constructed with exact 
analytic expressions for all components of the metric and fluxes. The background is compact 
and stable as a result of supersymmetry and a balance between fluxes and torsion. The 
geometry is complex, non-Kahler, and conformally Calabi-Yau. The wrapped D7-branes are 
located at the IR boundary where there are five internal angular directions. The induced 
fluxes warp space and result in four-dimensional spacetime at the UV boundary located at a 
finite distance from the IR boundary. The IR and UV boundaries are respectively located at 
r = r s and r = r s e 3 ^ on the radial direction R in^~S xR , where r s is the radius of the 
angular directions at the IR boundary and g s N is the 't Hooft coupling at the UV boundary. 
It was shown in [8] that the supergravity background reproduces the renormalization group 
flow and pattern of chiral symmetry breaking of the gauge theory and in [9] that it produces 
linear confinement of quarks within consistent relativistic quantum theory of ten-dimensional 
superstrings for the first time. In [10] mass spectrum of ++ glueballs in M = 1 SU(N) gauge 
theory in the large N limit was produced whose ratios happen to agree with available data 
from large N nonsupersymmetric lattice QCD to within 5%, which is consistent with linear 
confinement in both M = 1 supersymmetric and nonsupersymmetric SU(N) gauge theories. 

In this paper, we argue that the fact that the evolution of the universe cannot be ac- 
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counted by its visible matter and energy content might partly be due to forces external to 
it and propose a cosmological model of D3-brane universe on the supergravity background 
in [8] with the following assumptions: 

• A D3-brane universe of order Planck size filling four- dimensional spacetime R 1,3 spon- 
taneously appears in the UV region. 

• When the D3-brane collides with the wrapped D7-branes at the IR boundary, its kinetic 
energy is converted to matter and radiation confined to the D3-brane. 

The first assumption on spontaneous appearance of the model universe is similar to the case 
in the standard or other cosmological models where the origin or existence of a universe 
is always assumed. The second assumption is analogous to the process of reheating whose 
details are not well understood and is currently an assumption in inflationary cosmology 
models. 

The dynamics of the D3-brane until collision is governed by the Dirac-Born-Infeld (DBI) 
action that is determined by the background. The evolution of the D3-brane is analyzed 
starting with one located near the UV boundary. For examples of cosmological models and 
analyses which make use of the DBI action, see [11-21] for instance. The background we 
have here is compact and it does not involve AdS?, space since it is related to a confining 
gauge theory that is scale dependent. Compactness of the background geometry, availability 
of wrapped D7-branes at the IR boundary, and confinement of quarks that it produces 
provide a convenient setting for the beginning, reheating, and final stabilization of the model 
universe which then follows Friedmann evolution and it does not need to be glued to a larger 
background. Moreover, for N » 1 string loop corrections are negligible and the background 
is not measurably perturbed by the introduction of a single D3-brane. For small enough g s N, 
the curvature is small and a' corrections are negligible. 

It is found that the background geometry is such that the D3-brane flows toward the IR 
with accelerating expansion of its spatial volume followed by smooth transition to deceler- 
ating expansion. The cosmic evolution is dynamically dictated by the warped background 
geometry of the wrapped D7-branes with no additional branes or potentials introduced by 
hand. The total energy density of the D3-brane stays constant with the initial potential 
energy density getting converted to kinetic energy density as it flows toward the IR. Because 
the background preserves M = 1 supersymmetry in four dimensions, both the wrapped 
D7-branes and the probe D3-brane stay flat on R 13 . The magnitude of the scale factor of 
the three-dimensional volume produced is large enough to address the horizon and flatness 
problems provided that the D3-brane starts from close enough to the UV boundary. 
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The D3-brane reaches the IR boundary with velocity approaching its speed limit and col- 
lides with the wrapped D7-branes. We propose that the kinetic energy of the model universe 
is converted to matter and radiation by the collision followed by formation of baryons with 
strong interaction modeled by strings between quarks and the wrapped D7-branes serving 
as sources of color. The separation between quarks in baryons does not grow as the universe 
expands which gives a corresponding finite distance between the strings which end on quarks 
on the D3-brane and the D7-branes at the IR boundary, and strong interaction stabilizes the 
model universe at a finite distance from the wrapped D7-branes against gravitational force 
from the background. Friedmann evolution then takes over with a positive cosmological term 
coming from the remaining potential energy density of the model universe in the background 
which serves as a source of dark energy. The magnitude of dark energy density is smaller 
than the total energy density during the flow by a factor of the ratio of the scale factor when 
the model universe appears in the UV to the scale factor at the moment of collision at the 
IR boundary and stays constant while the matter-radiation density falls during Friedmann 
expansion, and it produces accelerating expansion during dark energy dominated phase. 

The inflationary scenario is used to generate sources of thermal fluctuations in the CMBR 
and seeds for large scale structure formation from quantum fluctuations which exit the Hub- 
ble radius early during the accelerating expansion phase, and the model addresses the inho- 
mogeneity problem. The power spectrum has red tilt consistent with observations. We show 
that slow-roll inflation can be used in this very early phase of cosmic evolution. Using results 
form slow-roll inflation and constraints from experimental cosmological observations on den- 
sity perturbations, spectral tilt, and dark energy together with measured value of Newton's 
gravitational constant in four dimensions, we write down constraints on the parameters of 
the model. 

We then consider one example of numerical values of the parameters which satisfy all the 
constraints except that the dark energy density is not small enough (but reduced from Planck 
scale in the right direction) and discuss their physical interpretations. Further investigation 
of the model, its parameters space, and available cosmological data is suggested considering 
the importance of addressing these problems simultaneously in the same setting which also 
provides a dual gravity theory for investigating the nonperturbative dynamics of Af = 1 
supersymmetric SU (N) gauge theory. We also discuss how the dynamics might be used to 
construct a cyclic universe. 
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2 Dynamics of D3-brane 

In this section, we write the DBI action for a probe D3-brane on the supergravity background 
of [8] and analyze its dynamics. First the ten-dimensional metric is 



dsl = r 2 ( cosh u g iiv dx ii dx v + sech u {dp 2 + dip 2 + dip 2 + dip 2 , + dip\ + dipl) ) , (2.1) 



x^ = (t,x) denotes the coordinates on four-dimensional spacetime with metric having 
(—,+,+, +) signature and r being time, (i/),ipi, ^2,^3,^4) are angular coordinates on the 
compact transverse internal space, and r is the radial coordinate. All components of x^ 
and r are measured in units of r s , the radius of the angular directions at the IR boundary. 1 
The coordinates in ten dimensions and the corresponding metric will be denoted by X M 
and Gmn- The dilaton $ is constant and the string coupling is given by = g s . The IR 
boundary is at r = 1 or p = and the UV boundary is at r = e 3 ^ or p = jj^, where g s N 
is the 't Hooft coupling at the UV boundary. The metric together with the fluxes on the 
background (all Fi, F3, if 3, and F§ fluxes of type IIB theory are turned on) solves all the 
supergravity equations of motion. See [8] for details on the supergravity background. 

2.1 Dirac-Born-Infeld action 

We want to study the dynamics of a probe D3-brane that fills four- dimensional spacetime 
and whose radial location on the background is represented by p = p(x M ). The metric given 
by (2.1) then becomes 



This notation avoids factors of r s floating around in our calculations and physical quantities are repre- 
sented by dimensionless variables in units of r s to appropriate power with c = ft = 1. We will put r s back 
and restore the dimensions to all variables, again in c = ft = 1 units, starting from the last paragraph of 
section 3.3 where we begin testing the model quantitatively with experimental data. 



where 




p = 3 In r, 



(2.2) 




(2.3) 
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The DBI action for the D3-brane is 



Sbbi = -~ / d^e-V-detGW, ( 2 - 4 ) 



where 

G^ u = r 2 cosh u(g^ + d^p d v p sech 2 u) (2.5) 

is the induced metric on the worldvolume and T3 is the tension of the D3-brane, measured in 
units of r~ 4 . Because there is no C4 potential in the background, there is no corresponding 
Chern-Simons term in (2.4). With (2.5), the action becomes 

Sdbi — ~ J d A x cosh 2 w \J — det((?^ + d^p d v p sech 2 u). (2.6) 

For the radial time evolution of the D3-brane excluding perturbations, which come in when 
we study inhomogeneities, p = p(r) and the action reduces to 

n 



Sdbi = / d x cosh u y 1 — p 2 sech u. (2.7) 

9s 



Extremizing (2.7), the equation of motion of the D3-brane is 

p = -{^Lf p cosh 4 u (2 - 3p 2 sech 2 u) . (2.8) 

271 

Rewriting the metric given by (2.3), now with p = p(r), as 

ds\ = r 2 coshw (1 — p 2 sech 2 u) {—dr 2 + a 2 dx^j 

+r 2 sechw (dip 2 + dy\ + dip 2 2 + dip 2 3 + dtpl) , (2.9) 



we see that the scale factor for the three-dimensional space is given by 



1 — p 2 sech u 



(2.10) 



Consider a D3-brane that is located near the UV boundary at p{r^) = — 5 at time 
Tj, where 5 « j 2 ^. We see from (2.8) and (2.2) that the D3-brane flows toward the IR 
boundary. Its initial acceleration is large provided that the initial location is close enough to 
the UV boundary, since coshu >> 1, and the evolution of the D3-brane is not sensitive to 
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its initial speed which we set to zero, p{ti) = 0. The scale factor given by (2.10) increases as 
the D3-brane flows toward the IR with the value of p getting close to its speed limit of — 1 
and sechw = 1 at the IR boundary. The magnitude of the scale factor when the D3-brane 
reaches the IR boundary can be made as large as desired by taking g s N 5 to be appropriately 
small. In the limit g s N 5 — > 0, the scale factor becomes infinitely large. The corresponding 
volume of the three-dimensional space grows by a factor of order ( ^g ) 3 - The D3-brane 
then collides with the stack of N wrapped D7-branes when it reaches the IR boundary at 



Numerical solutions to (2.8) are shown in figure 1 for the time evolution of p and a and 
their derivatives with p(rj) = j 2 ^ — 5, p(Tj) = 0, and the initial time Tj set to zero for 
g s N = 1 and 5 = 10~ 3 , with the value of the parameters chosen only for illustration of 
the features. We will write analytic expressions in section 2.2. Notice that the accelerating 
expansion takes place in a small fraction of the time interval during which the scale factor 
does not grow much. The speed of expansion then decreases while still staying large (>> 1) 
for nearly all of the flow for the parameters in the plot, and it is during this decelerating 
expansion phase that most of the growth of the scale factor occurs. The whole evolution is 
dynamically generated by the DBI action which is determined by the background; we have 
not introduced any additional terms to the equation of motion by hand. 

Thus the supergravity background produces dynamical accelerating expansion of the D3- 
brane followed by smooth transition to decelerating expansion with no additional ingredients 
introduced by hand, and the scale factor becomes larger as the initial location of the D3-brane 
gets closer to the UV boundary for a given g s N . 

Finally, the equation of motion of the D3-brane given by (2.8) can also be rewritten in 
different but equivalent forms such as 



p = 0. 



u = — 




2tt ' 



) 2 cosh 6 u — 5m 2 ) tanh u 



(2.11) 



which involves only u, and 




(2.12) 
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T 



Figure 1: Time evolution of p and a for the flow from near the UV boundary to the IR 
boundary and the corresponding p, a, p, and a for g s N = 1 and 5 = 10~ 3 . Some of the plots 
are made for shorter ranges of time to show the transition from accelerating to decelerating 
expansion closely. 
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2.2 Energy density 

The energy (or Hamiltonian) density for the D3-brane corresponding to the DBI action given 
by (2.7) is 

n= ?3 COs h^i = L cosh 2 M+( ^ )2cosh 4 M; (213) 

9s y/l - p2 S ech 2 u V 9s 
where p p is the conjugate momentum to p, 

p p =^ - t = . (2.14) 
y 1 — p 2 sech 2 m 

The potential and kinetic energy densities of the D3-brane are respectively given by 

V= — cosh 2 w (2.15) 

and 



T = — — 7^^^=^^^= cosh 2 -u = \ pi cosh 2 w + (— ) 2 cosh 4 -u cosh 2 -u. (2.16) 

9s y/l - p 2 sec h 2 u 9s V 9s 9s 

Using the energy density given by (2.13) and the equation of motion given by (2.8), we 
obtain 

f = 0. (2,7) 

Thus the total energy density of the probe D3-brane stays constant during the flow (and 
expansion) between the UV and IR boundaries. 

Consider a D3-brane that starts with piji) = at piji) = — 5 with 5 « j^. The 
energy density is %(rj) = cosh 2 u(xj) initially at time Tj and H(r c ) = y- — j= 1 ^ ^ when 
the D3-brane reaches the IR boundary at time r c . The energy density of the D3-brane stays 
constant with most of the energy density which initially at time is all potential converted 
to kinetic energy density at time r c , where T ~ p p . Because we take 5 « j 2 ^, we will from 
now on mostly do calculations to leading order in 9a ^ s and 

«(Ti) = --^ = tt(r), (2.18) 
9s 9sNb 

where H{r) is the energy density at any time between n and r c . Figure 2 shows the time 
evolution of the energy densities for g s N = 1 and 5 = 0.1, with the parameters chosen only 
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Figure 2: Total, kinetic, and potential energy densities of D3-brane measured in units of — 
for g s N = 1 and 5 = 0.1 with the first plot showing the time evolution until collision and 
the second one showing a shorter time range for the early phase. 



for illustration features including the remaining potential energy density when the D3-brane 
reaches the IR boundary. 

Analytic expressions for the scale factor a and the rate of change of p at any location p 
to leading order in are 

obtained using constancy of energy density, %{r) = "H(rj), 



a= — — — sech 2 w, p = — coshw \1 — (— cosh 2 u) 2 , (2-19) 

g s N6 V 7T 



1 1 ,g s N5 2 

- tanh up = — sinhu \ 1 — ( cosh u) 1 . (2.20) 

V 7T 



Notice that initially a(r») = 1 and p{ji) = as we had required and 
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When the D3-brane reaches the IR boundary, 



and the scale factor can be made large by taking g s N 5 to be appropriately small. The speed 
of the D3-brane when it collides with the wrapped D7-branes is such that p(r c ) -> -1 as 
g s N 5 — > 0. Because 5 << j 2 ^, p(r c ) f» —1. In fact, the D3-brane approaches the speed limit 
and p(r) ~ — 1 for most of the flow after the transition from accelerating to decelerating 
expansion. Notice that the scale factor at time of collision is related to the total energy 
density, 

o(r e ) = (2.23) 

Thus both the scale factor and the energy density increase with decreasing 5 as the initial 
position of the D3-brane gets closer to the UV boundary and become infinite when 5 = for 
a given g s N. We study the evolution of the D3-brane starting from 5^0, the background 
geometry is singular at the UV boundary where the curvature diverges because of warping 
of the transverse internal space to zero-size. 



3 Cosmological model 

In this section, we use the D3-brane dynamics studied in section 2 to construct a cosmological 
model and discuss its features. 



3.1 Expansion of model universe 

We assume that a model D3-brane universe of order Planck size spontaneously appears near 
the UV boundary at pfa) = — 5 at time r< with pfa) = 0, where 5 « j 2 ^. We have 
already studied the expansion in detail in section 2. Our convention for the scale factor is 
such that a(i~j) = 1. The model universe initially has potential energy and dynamically flows 
to the IR boundary while undergoing accelerating expansion followed by smooth transition 
to decelerating expansion and collides with the wrapped D7-branes at the IR boundary at 
time t c with speed p(r c ) = —yl — ( 2s ^) 2 and a scale factor of a(r c ) = j-jfg- The flow and 
expansion takes place with the total energy density staying constant, % = — JJv$- 

The kinetic energy of the D3-brane is converted to matter and radiation by the collision, 
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and using (2.16) and (2.18), 



n MR (r c ) = T(r c ) = ^ ( 1 - 1 ) = - - l ) ■ (3- 1 ) 

9s \y/l-p(r c ) 2 J 9s \9sN6 J 

Because the collision occurs instantaneously, the scale factor after collision equals the scale 
factor before collision. In order to obtain a scale factor that accommodates the observable 
size of the universe starting with order Planck size, we need -^ttt > 10 28 . 



3.2 Stabilization of model universe 

Now we use results in [9] to argue that the model universe bounces and stabilizes at a finite 
distance close to the IR boundary. When the kinetic energy of the model universe is converted 
to matter and radiation by the collision, quarks on the D3-brane begin forming bound states 
of baryons with strings in-between and the wrapped D7-branes serving as sources of color. 
Once baryons are formed, the size of the universe does not affect the size of the baryons; 
protons and neutrons continue to evolve in time with size of about a fermi. This is the 
scale of spatial separation L between quarks in baryons. For given parameters g s N and L, 
the profile of strings between quarks and their closest distance to the D7-branes at the IR 
boundary is finite. The model universe slightly bounces and stabilizes at a finite distance 
from the IR boundary at p = p e at time r e with a(r e ) « a(r c ) = j-jfg anc ^ strong interactions 
(or the tension in the strings) overcoming the gravitational force from the background. 

Let us briefly summarize the features we need from [9], with the focus here being a 
D3-brane located close to the IR boundary at p = p e << than at the UV boundary. 
See [9] for more details. Consider a quark-antiquark pair located at (x , x 2 , x 3 ) = (x, 0, 0) = 
(±7jr, 0, 0) on the D3-brane filling four-dimensional spacetime between the IR and UV bound- 
aries with a string in-between. The string worldsheet action for a static configuration at some 
time X° — > oo follows from the metric given by (2.1), and using Euclidean signature, 

S2 = *L i J da ° d(jl ^ detG MNd a X M d b X N = ^ J dx ^ cosh 2 u + , (3.2) 

where (a , a 1 ) parameterize the string worldsheet and we choose cr° = a; and a 1 = x 1 = x, 
p = p{x), and use the identification r s = yfot' in [8]. Minimizing the action (area bounded 
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by a rectangular Wilson loop) using the Euler-Lagrange equation, 



| = sech M cschn=^(i-(^) 2 p), (3.3) 



S 2 = 2 ( —(x - -sech 2 M )|^ = — L, 
\2n K 2 Jlx -°J 2tt ' 



(3.4) 



where we have taken p e « here (for the D3-brane near the IR boundary). Thus the 
potential energy of the quark-antiquark pair increases linearly with the spatial distance (or 
the string length) between them in four dimensions, E = For contrast, E = ^^p- when 
the D3-brane is at the UV boundary, which is the same as the energy we have here for 
L » 1. 

Solving (3.3) with boundary condition that x = ±4 at p = p e , we write the static 
configuration of the string for the region x > 0, it is symmetric about x = 0, as 



2n ^-i^V-pg)], (3-5) 



g a N\~ s po' 2 v 2tt 
where 

L 2tt p e lg s N 2 2 

= -t;M- ) - Ty^rr {Pe -Po ) (3-6) 
2 g s Jy po 2 27r 

and po measures how far the middle of the string is from the wrapped D7-branes. This is 
exactly the same result as in [9] for a D3-brane at the UV boundary which corresponds to 
p e = y^j. For given finite values of L and g s N, the value of ^ is finite too. The string 

ends on the quark and the antiquark with slope ^| p=Pe = (^-) 2 p e ) which has 

magnitude of -^h— for p e « we consider here and zero for p e = and 4^L-n = oo. 
The strings ending on quarks on the D3-brane then serve as anchors which support the 
model universe against gravitational force from the background, recall that most of the mass 
of the visible matter in the universe is in the strong force (or in the strings in this picture). 
Quarks in baryons move on the D3-brane universe with mean separation that does not change 
as the universe evolves. Schematic representation of a string in a quark-antiquark pair is 
shown in figure 3 based on the above results. Thus the model universe is stabilized against 
gravitational force from the background at a finite distance from the IR boundary by strong 
interaction which accounts for most of the visible mass of the universe. 
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Figure 3: Schematic description of the strings between quarks which stabilize the model 
D3-brane universe against gravitational force from the background of wrapped D7-branes. 
The D3-brane spontaneously appears in the UV region at p — j 2 ^ — 5 and dynamically flows 
toward the IR and collides with the D7-branes at the IR boundary at p = 0. The kinetic 
energy of the D3-brane is converted to matter and radiation by the collision and it stabilizes 
near the IR boundary at p = p e with the formation of baryons (a meson is shown here). 

3.3 The horizon, flatness, and inhomogeneity problems 

Now we discuss how the model addresses the horizon, flatness, and inhomogeneity problems. 
First let's briefly recall the nature of these problems in the hot big bang model and how 
the inflationary scenario addresses them. According to data from seven year Wilkinson 
Microwave Anisotropy Probe (WMAP) observations (combined with additional external data 
and fitting based on the standard cosmological model) [22], the age of the universe is 13.75 ± 
0.11 billion years, the current total energy density is Q tot = (1. 0023!$ 0054) ^c, the baryon 
density is Vl b = (0.0456 ± 0.0016) ft C) the dark energy density is fi A = (0.728!^) the 
dark matter density is Q^m = (0.227 ± 0.014) Q c , where Q c is the critical density which 
produces a flat universe in the Friedmann model, and the temperature power spectrum of 
the CMBR has a tilt parameterized by a scalar spectral index n s = 0.963 ± 0.012. In 
addition, the size of the observable universe is about 46 billion light years and the CMBR is 
homogeneous and isotropic with mean temperature of about 2.7 °K and relative variations 
of order 5 x 10~ 5 . 

In the hot big bang model the universe originated with sub-Planckian size and expanded 
to its current size. The horizon problem arises because if the current observable universe 
is evolved backward in time using the Friedmann equations for a radiation dominated case, 
one ends up with over 10 28 causally disconnected Planck size regions which makes it hard to 
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understand the observed homogeneity and isotropy. The flatness problem arises because the 
kinetic energy of the matter in the universe is related to the rate of expansion, and a flat 
universe with very near critical total energy density requires a highly fine-tuned initial veloc- 
ity to an accuracy of over 10 _54 %. The inhomogeneity problem is related to understanding 
the source of the small fluctuations of order 5 x 10~ 5 . 

The inflationary scenario addresses these problems by positing that the early universe 
underwent a period of accelerating expansion (inflation) which stretches its initially order 
Planck size by a factor of over 10 28 . The inflation arises due to a scalar field rolling on a 
potential energy density V(4>) that produces negative pressure and the sources of inhomo- 
geneities are quantum fluctuations which get stretched during inflation and early phase of 
Friedmann evolution. There are important requirements that the potential needs to satisfy 
in the standard slow-roll inflationary scenario in order to address the horizon, flatness, and 
inhomogeneity problems, see [23] for instance, 

M = |M 2 ^|«1, (3.8) 

1 1 v 3/2 

8 H = ^_ , „ — 5xl0~ 5 , (3.9) 



where M 4 is the Planck scale in four dimensions expressed in units of r" 1 = a'~ 1 ^ 2 here 
and a prime denotes differentiation with respect to 0. In standard inflationary models, 
the potential needs to be flat enough for a wide enough range of (f) m order to stretch the 
scale factor by over 10 28 and, at the same time, it needs to vary slowly enough at the time 
when quantum fluctuations exit the Hubble radius in order to produce the desired level of 
inhomogeneities. The observed spectral tilt in the CMBR spectrum puts the constraint 

n s - 1 = 2rj - 6e « -0.04. (3.10) 

In our case, first we show that in the region where accelerating expansion of the model 
universe occurs, we can use the methods of slow-roll inflation. Let us see how long the 
accelerating expansion lasts using the equation of motion given by (2.8), since p = at the 
transition from accelerating to decelerating expansion, which gives 

P( T eae) = ~ COsh U , (3.11) 
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where r eae is the time at the end of the accelerating expansion phase. We have < 
p 2 sech 2 u < 1 during this time. Therefore, we can make a Taylor expansion of the action to 
leading order in p 2 sech 2 u in this region and (2.7) reduces to 

S=f*x (1-P 2 -V), (3.12) 



2 2 



s 




where V is given by (2.15). Therefore, we can use the methods of slow-roll inflation during 
the accelerating expansion phase where we identify the scalar field 



(3.13) 



in units of r^ 1 (recall that T 3 and x M are dimensionless and measured in units of r~ 4 and 
r s respectively in our notation; we have (f> = \J~^ P r s with p = 31n(^) when expressed in 
terms of dimenionful variables). The potential given by (2.15) can rewritten in terms of 
as 

V=— ^ . (3.14) 

^i-i(^) 2 

First, with appropriate choice of parameters such that ^ > 10 28 , the scale factor 
grows by a factor of over 10 28 before a big bang collision followed by formation of baryons 
and stabilization near the IR boundary where a hot universe begins Friedmann evolution 
with a scale factor a(r c ) = j-jf$ Thus the model addresses the horizon problem. 

Moreover, because the background preserves M = 1 supersymmetry in four dimensions 
by a balance between fluxes and torsion, both the model D3-brane universe and the wrapped 
D7-branes stay flat along R 3 , the model universe expands from order Planck scale size by a 
large factor, and the production of matter and radiation due to collision of the flat branes 
also occurs homogeneously except for fluctuations. Thus the model addresses the flatness 
problem. In addition, the speed p(r c ) (and the kinetic energy and rate of expansion) at the 
time of collision is insensitive to the initial speed p(xj). That is because sechw(rj) « and 
the value of p(rj) 2 sech 2 ufa) « 1 insensitive to the value of p(rj); the kinetic energy of the 
model universe depends on the parameter 5 which describes how close its initial position is 
to the UV boundary for a given g s N. 

Next we focus on the generation of inhomogeneities that serve as sources of thermal 
fluctuations in the CMBR and seeds for large scale structure formation. For 93 j v<5 << 1, the 
acceleration at the beginning is given by (2.21). With (2.21) and (3.11), we can estimate 
the order of magnitude for the duration of the flow in which the expansion is accelerating 
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with a > 0, r eae — n ~ y 2gS ^- S3 ■ We also estimate the order of magnitude of how much the 
scale factor grows during this time and the rate at which it is changing at r eae using (3.11) 
in (2.10) and (2.20), a ^ e ) ~ \/3 and a(r eae ) ~ \J ^gsNS^ ' Therefore, the scale factor does 
not grow very much during the phase of accelerating expansion, its enormous growth occurs 
in the decelerating expansion phase where a starts with a large value for appropriately small 
g s N 5 3 and stays large for much of the flow until the model universe gets close to the IR 
boundary. This is different from standard inflationary scenarios where the enormous growth 
of the scale factor is produced during exponentially accelerating expansion. 

Quantum fluctuations that serve as sources of thermal fluctuations in the CMBR and 
seeds for large scale structure formation are produced during the accelerating expansion 
phase. In this phase the Hubble radius H~ x decreases while the scale factor a stays nearly 
constant. Therefore, quantum fluctuations of order Planck size exit the Hubble radius at 
some time r, where r, < r < r eae when their physical wavelength \ p h ys ~ H~ l . In order for 
the fluctuations to get stretched and grow with the scale factor, they need to stay outside 
the Hubble radius. This occurs if H~ x = f < \hys ~ Because both the fluctuations and 
the model universe begin with order Planck size, we need - < 1 which using (2.20) becomes 
tanhix|p| < Notice that tanI f M starts with infinitely large value at time Tj and falls during 
the accelerating expansion phase as |p| increases, and because tanhu w 1 in the UV, we 
require |p(r eae )| > 5 so that fluctuation exit the Hubble radius. Furthermore, we want 
5 << 1 (and 5 < 1 is good enough) so that the fluctuations stay outside the Hubble radius 
during most of the expansion as the model universe gets close to the IR boundary where 
\p\ ~ 1, a — >■ and tanhw — > 0. In this case, the fluctuations stay outside the Hubble radius 
for almost all of the decelerating expansion phase, since d >> 1 during most of the flow until 
collision. Close to the IR boundary, tanh u — > and the fluctuations enter the Hubble radius. 
The collision converts the kinetic energy of the D3-brane due to its motion along the radial 
direction to matter and radiation moving in four dimensions confined to the model universe 
which produces a discontinuity in a while a stays continuous, and the fluctuation exit the 
Hubble radius again. Furthermore, because fluctuations with larger \ p h ys exit the Hubble 
radius earlier during the accelerating expansion phase and have larger amplitude when they 
reenter, the power spectrum has red tilt consistent with observations. A schematic plot 
which shows how the Hubble radius H~ 1 and the physical wavelength of fluctuations \ p h ys 
evolve in time is shown in figure 4. Thus the model addresses the inhomogeneity problem. 

Because the background geometry is compact, M4 can be expressed in terms of the volume 
Vq from the extra space and we have a finite Newton's constant in four dimensions. Using 
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Figure 4: Schematic description of how fluctuations exit and enter the Hubble radius. Quan- 
tum fluctuations with physical wavelength \ p h ys of order Planck size cross the Hubble radius 
H^ 1 at time r ex close to r, during the accelerating expansion phase as the Hubble radius 
decreases. The Hubble radius increases as the model universe approaches the IR boundary 
and the fluctuations enter the curvature at time r, ni for a brief time during which the scale 
factor barely changes just before collision (big bang). The kinetic energy of the D3-brane is 
converted to matter and radiation by the collision at time r c , the rate of expansion jumps 
from zero to a finite value while the physical wavelength of fluctuations \ p h ys ~ a is contin- 
uous, and the fluctuations exit the Hubble radius again. The Hubble radius then increases 
and the fluctuations reenter the curvature at a later time Tj n2 during Friedmann evolution. 

the metric given by (2.1), the ten-dimensional volume is 



2-jr 



V w = V, (27r) 5 r s 6 J"" ) Jl-(^ p ydp = V,V 6 , V 6 = A7c 5 g s N, (3.15) 

where V4 = J d A x \f—g is the volume of the 4D spacetime, with V4 and Vq measured respec- 
tively in units of r 4 and r®, and we use the mapping a 1 = r 2 s in [8]. We then have 

2 2V 6 g s N 
M 4 = fn \7 2 = Ta 2 2 ' 3 ' 16 

with M 4 measured in units of r" 1 . Thus the background gives a finite value of Newton's 
gravitational constant in four dimensions which reproduces the experimentally measured value 
with appropriate choice of parameters. A smaller value of g s N ('t Hooft coupling at the UV 
boundary) corresponds to smaller curvature in which the supergravity description is valid 
on the whole region of interest away from the UV boundary. 
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We can obtain the total time the model universe takes from time r, to the beginning of 
Friedmann evolution using the equation of motion. Because the range of p it moves across 
is -jpy — 5 ~ and the speed p ~ — 1 for most of the flow to leading order in 9a ^ S , the 
order of magnitude of the total time is 

r c -r^^-. (3.17) 
9sN 

Using the relation between u and p given by (2.2) and the definition of given by (3.13), 

u' = 9 -^Mcosh*u, (3.18) 

where again a prime denotes differentiation with respect to <f), and with the potential in 
(2.15) we have 

9sN , 3 . 



V' = — - — \ — cosh wsinhw, (3.19) 



(9sN) 



2 



V" = - v ^ ' (2 cosh 2u - 1) cosh 4 u , (3.20) 
27T 

£ = - 9 —M ^h2«, = (2cosh2« - l)coshV (3.21) 

V 2tt V T 3 'V 2tt 2 T 3 v ; ' v 1 

where the minus sign is because the model universe is flowing toward the IR along which 
the potential decreases. 

The valid parameters space satisfying the constraints given by (3.7)-(3.10), (3.16), and 
the level of expansion can be investigated numerically. Our interest here is to find analytic 
expressions to leading order in 2^ in the region the model universe undergoes accelerating 
expansion and quantum fluctuations exit the Hubble radius. This occurs near p = j 2 ^ — 5 
where, to leading order in 9 'S fS , 

g s g s N6 V g s g s N6 2 g s Nd s 



(3.23) 



V fg^l _g^2_ 

V V T 3 5 ' V T 3 5 2 ' 

With (3.22) and (3.23), the constraints (3.7)-(3.10) together with (3.16) become, putting 
r s back which we do from now on with containing r s and the metric (2.1) rewritten as 
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in [8] , ds\ Q = cosh u g flu dx' Jl dx 1 ' + r 2 s sech u {dp 2 + dip 2 + dip 2 + dip>\ + d<p\ + dip 

M\ g 



2 ) 



1 ' s «0, (3.24) 
4e«0, (3.25) 



~ 2 5 2 T 3 r 2 

2M 2 g s 
5 2 T 3 r 2 ' 



1 ' 5 ^~5x 10- 5 , (3.26) 



\/757rM 4 3 y g s N g 

Ml = -^ty ■ (3.27) 



The model has the following parameters: g s N, 5, T3, r s , and iV (or g s ). We will consider an 
example of numerical values of these parameters and test the model in section 3.6. 



3.4 Dark energy and Friedmann evolution 

Once the model universe stabilizes, it follows Friedmann evolution with a cosmological con- 
stant term coming from the remaining potential energy density of the D3-brane in the back- 
ground which we interpret as dark energy. When it reaches the IR boundary almost all of 
the initial potential energy density is converted to kinetic energy density and the remaining 
potential energy density can be read off from the Hamiltonian given by (2.15) with p « 



2tt 



A = ^ = 44^)- (3-28) 
9s a(r c ) 

Thus the potential energy density of the model universe gives dark energy density that is 
smaller than the initial potential (or total) energy density by a ratio of the scale factor at the 
beginning near the UV boundary to that at the time of collision at the IR boundary. 
We then use this in the Friedmann equations for a flat universe, 

H ' = 3W {CM + A) - (3 ' 29) 

? = -J^( ej , + 8p-JA), (3.30) 

where H = - is the Hubble constant, Em is matter density, and p is pressure. 

For time r > r e fa r c , the matter-radiation density falls with increasing spatial volume 
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of the model universe while the dark energy density stays constant. From (3.1), (3.28), and 
(2.22) we see that 

= ^5~ 1 = a(Tc) - 1 (3 - 31) 

Because the total amount of matter-radiation for time r > r c should stay constant and 
a(r) » 1 (in our notation) during Friedmann evolution, the ratio of the matter-radiation 
energy density now at time r n to that at the beginning of Friedmann evolution is 



UMR{r n ) f a(r c ) x " 



Hmr{t c ) \a(r n ) 



(3.32) 



On the other hand, cosmological observations show that the current ratio ^ M ^( T ") | ; where 
we have included cold dark matter in %mr- Thus the scale factor today is related to the 
scale factor at the beginning of Friedmann evolution, 

^Ktf /3aW M?Ffef • (3 - 33) 

According to the picture here, a constant dark energy density existed since collision 
(big bang), and its contribution to the evolution of the universe becomes proportionally 
more significant as the universe expands and the matter-radiation density become smaller. 
Thus the remaining potential energy density serves as a source of dark energy that produces 
accelerating expansion of the model universe during dark energy dominated phase. 



3.5 Summary of constraints from observations 

First let us recall the parameters of the model: g s N, 5, T3, r s , and N. The 't Hooft coupling 
at the UV boundary is g s N and its magnitude determines the range of p, < p < j 2 ^. A 
smaller value of g s N corresponds to a larger range and a smaller curvature. The radius of 
the internal space at the IR boundary is r s . N is the number of wrapped D7-branes at the 
IR boundary. The model universe spontaneously appears in the UV region at p = — 5 
at time Tj with 5 « j 2 ^. The initial energy density of the D3-brane is all potential and 
%{Ti) = y YWsi an d its magnitude is larger than the dark energy density A = ^ during 
Friedmann evolution by a factor of -^rn- 

J g s N8 

Let us summarize the constraints that follow from the desired level of expansion needed 
to address the horizon and flatness problems for a universe that begins with Planck size, 
cosmological data on inhomogeneities, the measured value of Planck scale (or Newton's 
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constant) in four dimensions, and observed level of current accelerating expansion, 



-jH = ~frri ~ 10 28 , (3.34) 
aiji) g s N 6 

Ml g s 1 



2 5 2 T 3 r 2 s 4 



V\ , (3-35) 



1 / 5 T-?r s c 

- 5 x 10" 5 , (3.36) 



Vl^Ml V 9sN g s 



J/i : \/ 1/2^2 2 = 1 « ^ 35 = L2 x 10 19 GW, (3.37) 
" lQn 2 g 2 r 2 1.6 x 10~ 35 m 

A= — ~ 1(T 120 M 4 4 . (3.38) 
The constraints on the values of e and rj gives the spectral tilt in the CMBR, 

n s - 1 = 2r] - 6e = -14 e w -0.04. (3.39) 

In addition, requiring that the fluctuations exit the Hubble radius during most of the expan- 
sion as shown in figure 4 requires 

5 < 1 . (3.40) 



3.6 An example of numerical values of parameters 

The model has already been constructed in previous sections. Here we discuss one example 
of numerical values of the parameters and test the model. 
Consider the following choice of parameters, 

g s N = 7.2 x 10~ 28 , T 3 = 5.9 x 10 7 r7 4 , M 4 = 3.0 x 10 19 r' 1 , 

S = 4.3 x 10- 1 , N = 10 7 , r s = 4.7 x KT 16 m. (3.41) 

First, notice that r s (the radius of the internal space at the IR boundary) is about the size of 
baryons which is a reasonable value for identifying r s with the nonperturbative scale of the 
gauge theory whereby strong interaction stabilizes the model universe near the IR boundary 
against gravitational force from the background. The volume of the extra space seen by a 
four-dimensional observer and given in (3.15) is Vq = A% 5 r & s g s N = (4.6 x 10~ 20 m) 6 . We 
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also have the measured numerical value of Planck scale (or Newton's gravitational constant) 
in four-dimensions M 4 = J ' a = 3.0 x 10 19 rr 1 = = - * , 5 = 1.2 x 10 19 GeU and 

* V loir^r-gg^ a 1.6x10 *m 

a desirable level of expansion 4^4 = T r x = 10 28 . It takes a time of about -^j— = 3.8 

1 a(Ti) g s N 6 g s N c 

hours for the model universe to flow from near the UV boundary to collision at the IR 
boundary which is taken as a big bang and the beginning of a hot universe with matter and 
radiation. The slow-roll constraints on the potentials given by (3.24)-(3.26) are satisfied, 
e = gjij = 2.9x 1(T 3 « 1 and = ^ ^ = 4e = 1.1 x 1(T 2 « 1. For the spectral 
tilt in the CMBR, n s — l = 2r] — 6e = —0.04 and for the density perturbation, 5h = 5 x 10~ 5 . 
Moreover, because 5 = 0.43 < 1, quantum fluctuations which exit the Hubble radius during 
the accelerating expansion phase stay outside and grow with the scale factor as schematically 
shown in figure 4. Thus the above set of parameters produces the desired level of expansion 
for a universe that begins with Planck size in addition to the observed experimental values 
of scalar density perturbation, spectral tilt, and Newton's gravitational constant in four 
dimensions. 

The dark energy density during Friedmann evolution is "Hde = ~" — 8.1 x 10 41 r~ 4 = 
1.1 x 10~ 36 M 4 . The total energy density, which stays constant during the flow until the 
beginning of Friedmann evolution and equals the potential energy density at time Tj, is 
U(t) = J ^ = 1.1 x 10~ 8 M 4 4 = (0.01 M 4 ) 4 for n < r c . Thus, for the numbers we have 
here, the model universe flows with order Planck total energy density which is initially all 
potential and almost all of it becomes kinetic at the IR boundary. The kinetic energy is 
then converted to matter and radiation by the collision. The remaining potential energy 
density becomes dark energy density that is smaller than ~H(t) by a factor of 10 -28 ; i.e., A = 
1.1 x 10 -36 M|, which is reduced in the right direction but not small enough to accommodate 
the observed value of order 10~ 120 M\. Recall that "Hde = A = ^\'H(t) and a smaller value 
of dark energy density can be obtained by taking smaller g s N 5. One needs ^|^y ~ 10 120 
with H{t) of order Planck energy density in order to address the dark energy problem, and 
this needs to be done while the other constraints on the desired level of expansion, density 
perturbations, and spectral tilt are satisfied. 

Finally, because g s N (the 't Hooft coupling at the UV boundary) is small, the curvature 
throughout the region in which the model universe flows is small and a' corrections are 
negligible. The string coupling at the UV boundary is g s = ^j^- = 7.2 x 10~ 35 , which is 
small. Neither the gauge theory nor the gravity theory predicts the absolute magnitude of 
the string or gauge coupling (the renormalization group flow relates the coupling at one scale 
to the coupling at another) and we have picked a value that is suitable for our discussion here. 
The string we have here is that of strong interactions with tension of = 5 where r s is 
given above. Furthermore the D3-brane tension T 3 = ( ^ r ) 4 , which has the same order of 



24 



mass scale. When the D3-brane appears in the supergravity background, its energy density 
becomes order Planck scale, as we saw above, due to gravitational interactions. The four- 
dimensional gauge coupling runs and both the gauge and 't Hooft couplings of the gauge 
theory living on the wrapped D7-branes become 0(1) and larger close enough to the IR 
boundary as r — > r s or p — > 0, see [8] for details. Moreover, because iV >> 1 both string 
loop corrections and backreaction from introducing a single D3-brane are negligible. We 
would also like to point out that lattice QCD computations for physical observables such as 
mass ratios of glueballs in the iV — > oo limit are not significantly different from the case of 
N = 3, see [24] for instance. 

4 Conclusions 

The fact that the evolution of the universe cannot be accounted by its visible matter and 
energy content might partly be due to forces external to it and the supergravity background 
argued to be dual to pure Af = 1 SU(N) gauge theory in [8-10] produces accelerating 
expansion followed by smooth transition to decelerating expansion on a D3-brane model 
universe near the UV boundary as it flows toward the IR with the dynamics dictated by 
the background geometry through the DBI action without introduction of additional terms 
by hand. The model presented addresses the horizon, flatness, and inhomogeneity problems 
and provides a possible explanation to dark energy. 

The model universe flows to the IR boundary with constant total energy density with 
almost all of the initial potential energy density becoming kinetic energy density. The kinetic 
energy is then converted to matter and radiation by collision (big bang) with the wrapped 
D7-branes at the IR boundary. This conversion of kinetic energy to matter and radiation 
is one of the assumptions we have made in constructing the model, and it is analogous to 
the process of reheating which is similarly an assumption in standard inflationary models. 
The remaining potential energy density gives a cosmological constant which is interpreted 
as dark energy density. The strong interaction which accounts for most of the visible mass 
of the universe is modeled by strings between quarks in baryons (whose size does not change 
as the universe expands) which stabilizes the model universe at a finite distance from the IR 
boundary against gravitational force from the background. The model universe then follows 
Friedmann evolution. Because the dark energy density stays constant while the matter- 
radiation density falls during Friedmann expansion, the contribution of dark energy to the 
cosmic evolution increases with increasing time and produces accelerating expansion in dark 
energy dominated phase. The ratio of dark energy density during Friedmann evolution to 
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the initial energy density when the model universe appears near the UV boundary (or during 
the flow) is equal to the ratio of the scale factor initially to the scale factor at the beginning 
of Friedmann evolution. 

Cosmological observations from the desired level of expansion for a universe that begins 
with order Planck size and experimental values of density perturbations, spectral tilt, dark 
energy, and Newton's gravitational constant in four dimensions put constraints on the pa- 
rameters of the model. We have written down the constraints using results from slow-roll 
inflation and discussed one example of numerical values of parameters which accommodates 
all the constraints except that the dark energy density is not small enough (but reduced from 
Planck scale in the right direction). Smaller value of dark energy density can be obtained 
by taking smaller g s N 5. Considering the importance of addressing the problems discussed 
here simultaneously in the same setting, in addition to providing a dual theory for studying 
the nonperturbative dynamics of M = 1 supersymmetric SU(N) gauge theory, it will be 
interesting to investigate the model, its parameters space, and available cosmological data 
further. 

Because the total energy of a Planck size model universe is comparatively much smaller 
than the energy content of a large number of wrapped D7-branes, it is plausible that the 
model universe might originate from quantum fluctuations of the background. The spon- 
taneous appearance of a D3-brane in the background is the other assumption we made in 
constructing the model; in the standard or other cosmological models, the origin or existence 
of a universe is always an assumption. 

The dynamics discussed in section 2 might also be useful to investigate possibilities for 
cyclic universe scenarios such as implementing some features of the ekpyrotic scenario [25] 
and its cyclic extension without requiring a singular beginning [26]. The size between the 
IR and UV boundaries and the time interval between collisions can be made large by taking 
g s N to be appropriately small. If only part of the kinetic energy of the D3-brane is converted 
to matter and radiation by the collision, then the D3-brane would continue bouncing back 
and forth (returning from a point which gets further away from the UV boundary after each 
bounce). See appendix A for more. 

Finally, the results in this paper together with [8-10] lead us to a picture of D3-brane 
model universe with matter and radiation on it, the wrapped D7-branes serving as sources 
of color to strong interactions, and dark energy interpreted as potential energy of the model 
universe in the background. Thus the supergravity background can be used to investigate 
both the nonperturbative dynamics of strong interactions and cosmology simultaneously in 
the same setting. 
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A Cyclic universe 



Here we discuss how the dynamics discussed in section 2 might also be useful to construct 
a dynamical cyclic universe. Consider a D3-brane that is located near the UV boundary 
at p(Tj) = j 2 ^ — 5 at time Tj with p(rj) = 0, where 5 << Tj does not need to be a 
beginning. The D3-brane flows toward the IR with the scale factor of the three-dimensional 
space given by (2.10) increasing to a(r c ) = j-jfg a(Tj) at time of collision at the IR boundary. 
Part of the potential energy density of the D3-brane is converted to kinetic energy density 
when the D3-brane reaches the IR boundary. Considering a case in which only part of the 
kinetic energy is converted to matter and radiation by the collision, the D3-brane bounces 
back with a reduced radial speed and comes to a stop at a further point from the UV 
boundary than its location at time n and then begins to flow back to the IR again. The 
cycle repeats. Thus the background provides a setting for cyclic universe scenarios. This 
might, for instance, be useful to investigate possibilities for implementing some features of 
the ekpyrotic model [25] and its cyclic extension [26]. Here, the D3-brane bouncing back 
and forth might be taken as the universe and its cosmic evolution is dynamically dictated by 
the warped background geometry of the wrapped D7-branes. The size between the IR and 
UV boundaries and the time interval between collisions can be made large by taking g s N 
to be appropriately small. The warping of space is such that the volume of the extra space 
seen by a four-dimensional observer on the D3-brane and the curvature of the background 
decease with increasing distance between the two boundaries (or decreasing value of g s N). 
The D3-brane does not need to be of Planck size and can be large when it is near the UV 
boundary or at time Tj. The four- dimensional spacetime here is smooth throughout with 
the singularity at the UV boundary arising from warping of the transverse internal space to 
zero-size. Plots for the time evolution of p and a for a case of elastic collision are shown in 
figure 5 for g s N = 10~ 3 and 5=1 with the time Tj set to zero and the parameters chosen 
only for illustration; the D3-brane expands as it flows toward the IR and contracts as it 
bounces back to the UV. 
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